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ABSTRACT 

The anisotropic galaxy 2-point correlation function (2PCF) allows measurement of the 
growth of large-scale structures from the effect of peculiar velocities on the clustering pat¬ 
tern. We present new measurements of the auto- and cross- correlation function multipoles 
of 69,180 WiggleZ and 46,380 BOSS-CMASS galaxies sharing an overlapping volume of 
~ 0.2 (/i -1 Gpc) 3 . Analysing the redshift-space distortions (RSD) of galaxy 2-point statis¬ 
tics for these two galaxy tracers, we test for systematic errors in the modelling depending 
on galaxy type and investigate potential improvements in cosmological constraints. We 
build a large number of mock galaxy catalogs to examine the limits of different RSD 
models in terms of fitting scales and galaxy type, and to study the covariance of the 
measurements when performing joint fits. For the galaxy data, fitting the monopole and 
quadrupole of the WiggleZ 2PCF on scales 24 < s < 80 h ~ x Mpc produces a measurement 
of the normalised growth rate f(J%(z = 0.54) = 0.409 4= 0.059, whereas for the CMASS 
galaxies we found a consistent constraint of fcr%(z — 0.54) = 0.466 ±0.074, When combin¬ 
ing the measurements, accounting for the correlation between the two surveys, we obtain 
fag(z = 0.54) = 0.413 d= 0.054, in agreement with the ACDM-GR model of structure 
growth and with other survey measurements. 

Key words: cosmology - large scale structure 


1 INTRODUCTION 

The evolution of the spatial distribution of galaxies on large 
scales is deeply influenced by the physics of gravitational at¬ 
traction, cosmic expansion, and the conditions of the early 
Universe, and constitutes an important probe and discrim¬ 
inator of cosmological models. Spectroscopic galaxy surveys 
map this distribution using the distance-redshift relation, but 
due to peculiar velocities induced by the gravitational field, 
these maps also contain ‘redshift-space distortions’ (RSD) of 
the spatial positions of the galaxies, which modify the true 
(i.e. real space) pattern of the spatial clustering of galaxies. 
Kaiser (1987) showed that on large scales the peculiar veloc¬ 
ity field v (in dimensionless units of the Hubble velocity) is 
related to the matter overdensity 5m as V • v = — f5 m , where 
the proportionality parameter f(z) is called the linear growth 
rate of structure. Modelling the redshift-space clustering, in 
consequence, allows us to constrain cosmological parameters 

* E-mail: fmarin@astro.swin.edu.au (FAM) 


through estimations of f(z). Using Kaiser’s findings, the pi¬ 
oneering works in the 2dFGRS survey (Peacock et al. 2001; 
Hawkins et al. 2003) in the local z ~ 0.1 Universe, measured 
the redshift-space two-point clustering of galaxies, which re¬ 
sulted in a confirmation of the concordance ACDM model at 
present times. With the advent of galaxy surveys at higher 
redshifts, we can now trace the history of f(z) and obtain 
constraints on cosmological models and the nature of dark 
energy (Linder & Cahn 2007). 

However, important challenges must be addressed before 
we can use this tool effectively. On the observational side, the 
most important factors limiting the statistical precision of the 
clustering measurements obtained from different surveys are 
the limited volume that the surveys can map, due to the sam¬ 
ple variance from fluctuations in the clustering on different 
regions of the universe, and the discreteness of the galaxy field 
known as shot noise (e.g. Kaiser 1986; White et al. 2009) 

In addition, large-scale structures are subject to a vari¬ 
ety of systematic non-linear effects which affect our capac¬ 
ity to model the signal, particularly on small scales. First, 
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we have non-linear growth of structure, such that even on 
large scales, the Kaiser relations are insufficient to account for 
the measured clustering in galaxy data and simulations. Pea¬ 
cock (1992), and more recently Scoccimarro (2004); Taruya 
et al. (2010); Seljak k McDonald (2011); Wang et al. (2013), 
among others, have improved the basic ‘Kaiser’ model by in¬ 
cluding various non-linear effects in the matter clustering. Sec¬ 
ond, there is scale-dependent complexity in how galaxies trace 
haloes and cross-correlate to matter, known as galaxy bias. 
Third, galaxies possess non-linear pairwise velocities on small 
scales. The latest attempts to use the 2-point clustering pat¬ 
tern to model RSD have taken these and other effects into ac¬ 
count (e.g. Reid et al. 2012; Beutler et al. 2012; Sanchez et al. 
2013; de la Torre et al. 2013; Contreras et al. 2013; Beutler 
et al. 2014, as recent examples), allowing us to confront pre¬ 
dictions from different cosmological models. 

Although in linear theory all galaxies respond as test par¬ 
ticles to the gravitational field, in detail the non-linear system- 
atics depend on tracers themselves. This is because galaxy for¬ 
mation is affected by many non-linear processes such as small- 
scale dynamics of halo formation, environment, and complex 
baryonic processes determining the luminosity and colour at 
a given time, which are the main observables when select¬ 
ing galaxies for a large-scale survey. Therefore, the analysis 
and modelling of two overlapping tracers makes it possible 
to constrain details of the clustering and formation of the 
galaxy tracers themselves. Previous work has focused on cross- 
correlating a tracer with well known properties with a sec¬ 
ond tracer we wish to study (e.g. Martinez et al. 1999; Chen 
2009; Mountrichas et al. 2009; Font-Ribera et al. 2013). In our 
current study we approach this in a cosmological context, in 
which a comparison of results using different tracers in the 
same volume tests for systematic errors in modelling of bias 
and redshift-space distortions. 

In this work we present measurements and analysis of 
RSD using galaxies from the WiggleZ Dark Energy Survey 
(Drinkwater et al. 2010) and the CMASS galaxy sample from 
the Baryon Oscillation Spectroscopic Survey (BOSS, Eisen- 
stein et al. 2011). At a redshift of z ~ 0.6, the WiggleZ team 
targeted Emission Line Galaxies hosted in low-to-intermediate 
mass halos, which have low bias (frwiggiez ~ 1, see Blake et al. 
2011b; Contreras et al. 2013; Marin et al. 2013), while the 
CMASS sample consists of luminous, mostly red galaxies with 
&CMASS ~ 2 (Reid et al. 2012; Tojeiro et al. 2012; Chuang 
et al. 2013; Kazin et al. 2013) with similar number density 
n ^ 2 — 3 x 10 -4 (/i -1 Mpc) -3 . With an overlap volume of ap¬ 
proximately 0.2 (h -1 Mpc) -3 , this is, to date, the largest vol¬ 
ume overlapping sample between two different galaxy redshift 
surveys. We measure the redshift-space auto- and cross- corre¬ 
lation functions of these galaxies and explore the constraints 
on the cosmic growth rate using the two tracers. Our work is 
supported by a large suite of mock catalogs, which we gen¬ 
erated by performing abbreviated N-body methods (COLA, 
Tassev et al. 2013) to model potential systematics coming 
from observational issues, test different RSD models and their 
regime of validity, and determine covariances. 

A potential advantage of a multi-tracer analysis was de¬ 
scribed by McDonald k Seljak (2009), who noted that the 
correlations in an overlapping volume, if the number density 
of the tracers is large, can be used to reduce the sample vari¬ 
ance error and improve the measurements of the growth rate. 
After this initial work, different applications of the multitracer 
method have been explored by various authors, using different 


observables such as photometric redshift surveys, weak lens- 
ing, gravitational redshifts, signatures of first stars and con¬ 
straints on primordial non-gaussianity and modified gravity 
(Seljak 2009; Bernstein k Cai 2011; Gaztanaga et al. 2012; 
Asorey et al. 2013; Croft 2013; Yoo k Seljak 2013; Lombriser 
et al. 2013). Blake et al. (2013) applied this method to the 
GAMA survey, producing modest gains from the multitracer 
method, up to 20% in the constraints of / at two different 
epochs, z = 0.18 and z — 0.38. Ross et al. (2014) measured 
the clustering of BOSS galaxies as a function of their colour 
and did not detect significant differences in distance scale or 
structure growth measurements. Although the datasets used 
in our study are too sparse to expect large improvement, we 
include this effect by computing the full covariance of the mea¬ 
surements using our mock galaxy catalogs. 

We present in section §2 the surveys used in our study. 
In §3 we present the methods and results of the auto- and 
cross- correlation between tracers. In §4 we show models of 
the RSD and constraints in the model parameters and the 
growth rate at z — 0.54. Finally in §5 we summarize our 
results and conclude. This is the second work of a series of 
papers analysing clustering in the BOSS-WiggleZ overlap re¬ 
gion. Paper I (Beutler et al., 2015) focuses on the analysis of 
the Baryonic Acoustic Oscillation signal of these two tracers 
in the common volume. 

For clarity we will use the name ‘CMASS-BW’ and 
‘WiggleZ-BW’ for the CMASS and WiggleZ samples limited 
to the overlap region between the two surveys. We assume a 
fiducial flat ACDM cosmological model as defined in Komatsu 
et al. (2009), where the matter density is Q m = 0.273, baryon 
density of Ob = 0.045, a spectral index of n s — 0.963, an 
r.m.s. of density fluctuations averaged in spheres of radii at 8 
h _1 Mpc of <7 8 — 0.81 and h — 0.71. The Hubble rate at red¬ 
shift z =0 is Ho= 100/i km s -1 Mpc -1 is adopted to convert 
redshifts to distances, which are measured in /i _1 Mpc. 


2 DATA k MOCK CATALOGS 
2.1 The WiggleZ survey 

The WiggleZ Dark Energy Survey (Drinkwater et al. 2010) is 
a large-scale galaxy redshift survey performed over 276 nights 
with the AAOmega spectrograph (Sharp et al. 2006) on the 
3.9m Anglo-Australian Telescope. With a area coverage of 816 
deg 2 , this survey has mapped 207,000 bright emission-line 
galaxies over a redshift range 0.1 < z < 1.0. Target galax¬ 
ies in six different regions were chosen using UV photometric 
data from the GALEX survey (Martin et al. 2005) matched 
with optical photometry from the Sloan Digital Sky Survey 
(SDSS DR4, Adelman-McCarthy et al. 2006) and from the 
Red-Sequence Cluster Survey 2 (RCS2, Gilbank et al. 2011). 
The selection criteria consisted of applying magnitude and 
colour cuts (Drinkwater et al. 2010) in order to select star¬ 
forming galaxies with bright emission lines with a redshift 
distribution centered around z ~ 0.6. The selected galaxies 
were observed in 1-hour exposures using the AAOmega spec¬ 
trograph, and their redshifts were estimated from strong emis¬ 
sion lines. The number density of WiggleZ galaxies averages 
~ 3 x 10 -4 (/i -1 Mpc) -3 at z = 0.6. 
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Figure 1 . Sky coverage of BOSS-CMASS DR11 (red) and WiggleZ (blue). The left panel shows the northern part of the surveys, while the 
right panel shows the southern sky coverage. Five of the six WiggleZ regions are covered by CMASS, with a total of 69,180 WiggleZ galaxies 
and 46,380 CMASS galaxies in the overlap volume. 


2.2 The CMASS Sample 

The Baryon Oscillation Spectroscopic Survey of the Sloan Dig¬ 
ital Sky Survey III (SDSS-III, Eisenstein et al. 2011; Daw¬ 
son et al. 2013), which is now complete, was designed to ob¬ 
tain spectra and redshifts for 1.35 million bright galaxies over 
a footprint ~ 10,000 deg 2 . These galaxies are selected from 
the SDSS-III imaging and have been observed together with 
160,000 quasars and 100,000 ancillary targets (Gunn et al. 
2006; Bolton et al. 2012; Smee et al. 2013). The CMASS sam¬ 
ple is composed of luminous, mostly red galaxies selected to 
probe large-scale structure at intermediate redshifts, achiev¬ 
ing a number density of ~ 3 x 10 -4 (/i -1 Mpc) -3 . The DR11 
catalog (Alam et al. 2015) includes 1,100,000 spectra out of 
which the CMASS sample contains ~ 550,000 galaxies in the 
redshift range 0.43 < z < 0.7. 

2.3 Overlap volumes 

We define the overlap regions between CMASS and WiggleZ 
using the random galaxy catalogs generated for each survey, 
gridding the sky into 0.1 deg 2 regions and selecting cells con¬ 
taining both CMASS and WiggleZ random points. As seen 
in Figure 1, five of the six WiggleZ regions have considerable 
overlap with CMASS galaxies, totalling 560 deg 2 and a vol¬ 
ume of 0.218 (/i -1 Gpc) 3 in the 0.43 < z < 0.7 range. This 
results in an overlap sample of 69,180 WiggleZ galaxies and 
46,380 CMASS galaxies. 

Figure 2 shows the redshift distribution of the two sam¬ 
ples in the different regions, which is similar in the range 
0.5 < z < 0.6; outside that range the CMASS galaxy counts 
rapidly decline. To estimate how these differences will affect 
our results, we calculate the pair-weighted redshift, which con¬ 
sists in taking the average redshift of all pairs at a particular 
distance range. For the distance range s = 20 — 35 /i _1 Mpc, 
where the signal of the clustering signal is higher, WiggleZ-BW 
galaxies have a pair weighted redshift of £wiggieZ-Bw = 0.56 
whereas for CMASS-BW galaxies zcmass-bw = 0.53. For 
cross-pairs this redshift is z aV g,x — 0.54. These small differ¬ 
ences in redshift will not affect our findings given the mea¬ 
surement errors, therefore we generate cosmological models at 



redshift z redshift z redshift z 

Figure 2. Number of galaxies as a function of redshift for WiggleZ 
(blue) and CMASS (red) galaxies in the overlap regions, and for the 
overall overlap volume (bottom right panel). 

z a vg,x = 0.54 to compare with our WiggleZ-CMASS cluster¬ 
ing data. Table 1 presents details of the samples used. 

2.4 Simulations and mock catalogs 

We estimate the covariance of our measurements and test the 
regime of validity of our RSD models using mock galaxy cata¬ 
logs built from N-body simulations. The conventional methods 
to generate N-body simulations do not allow for the genera¬ 
tion of a large number of realisations of cosmological volumes 
at sufficient mass resolution to encompass the low-mass halos 
hosting WiggleZ galaxies, which are needed for constructing 
robust covariance matrices. For this reason we use an approx¬ 
imate, fast method to generate dark matter simulations based 
on the COmoving Lagrangian Acceleration method (COLA, 
Tassev, Zaldarriaga, & Eisenstein 2013). We have developed 
a parallel version of COLA (Koda et al., in preparation, used 
first in Kazin et al. 2014), where in each simulation contains 
1296 3 particles in a box of side 600/i _1 Mpc, which gives a par¬ 
ticle mass of 7.5 x 10 9 /i _ 1 M© , allowing resolution of low-biased 
halos with masses 10 12 /i _1 M©, found using friends-of-friends 
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Table 1. Overlapping samples analysed in this study. 


Region 

WiggleZ-BW 

N aal 

WiggleZ-BW 

Zavg 

CMASS-BW 

N qal 

CMASS-BW 

Zavg 

Cross-pairs 

Zavg 

SOI 

6620 

0.61 

5720 

0.53 

0.57 

N09 

13940 

0.56 

9360 

0.53 

0.54 

Nil 

15560 

0.55 

10580 

0.53 

0.54 

N15 

22740 

0.56 

14660 

0.54 

0.54 

S22 

10320 

0.55 

6060 

0.53 

0.54 

Total 

69180 

0.56 

46380 

0.53 

0.54 



Figure 3. Projected correlation function w p (r p ) for WiggleZ and 
CMASS galaxies in the overlap regions (symbols). Lines are the 
mean values of w p (r p ) for the COLA mock catalogs. 


algorithm with a linking length of 0.2 times the mean parti¬ 
cle separation. Each simulation requires 15 minutes with 216 
computation cores, including halo finding, which is much faster 
than a classical N-body simulation, but with similar precision 
on the relevant scales {k < 1 /iMpc -1 ). 

We generate a total of 2400 realisations (480 for each 
WiggleZ region) of a flat ACDM universe with WMAP5 cos¬ 
mological parameters (Komatsu et al., 2009), which defines 
our fiducial cosmology. Using the output at z = 0.6 we cre¬ 
ate WiggleZ-based (WiZcola) and CMASS-based (BOSScola) 
mock galaxy catalogs, from simple Halo Occupation Distri¬ 
bution models (Berlind & Weinberg 2002; Blake et al. 2008), 
such that the resulting projected correlation functions w p (r p ) 
match those of the observations, as seen in Figure 3. We then 
apply the relevant selection functions to the mock galaxies to 
match the survey geometry. Our simulations encode the joint 
covariance in the overlapping survey regions (Koda et al., in 
prep.). 


3 MEASUREMENTS 

3.1 Measuring Correlation Functions 

We estimate the redshift-space two-point correlation function 
£(s,/x) (2PCF) as a function of comoving separation s and 
the cosine of the angle of the distance vector with respect 
to the line of sight n — cos (0). We use the Landy &; Szalay 
(1993) estimator, counting pairs of objects in data and random 


catalogs: 

^ ^ ^ _ DD(s,n)-2DR(s,n) + RR{s,n) ^ 

RR(a,n) ’ 

where PP, DR and RR are respectively the weight- 
normalised data-data, data-random and random-random pairs 
with separation s and i± (with a given resolution As, A fi). 
For both the random and data catalogs we use the optimal 
(inverse-density) FKP weighting (Feldman et al. 1994): 


Wi(x) 


1 

1 + n(x)P 0 


( 2 ) 


where Po = 5000 (h -1 Mpc) 3 for WiggleZ-BW and Po — 20000 
(h _1 Mpc) 3 for CMASS-BW galaxies. For WiggleZ galaxies, 
angular incompleteness and radial selection are introduced in 
the random catalogs (Blake et al. 2010). A small fraction of 
galaxies contain errors in the redshift assignment, but this 
effect is absorbed into the fitted galaxy bias factor. CMASS 
galaxies, have additional weights applied to account for the an¬ 
gular incompleteness, fibre collisions, redshift failure and cor¬ 
relation between density of targets and density of stars (Ross 
et al. 2012). 

It is possible to model the 2PCF using the full information 
from £(s,/x), but that requires a large covariance matrix with 
the associated problems with its inversion. For this reason it 
is standard to compress this information in multipoles 




21 + 1 
2 



£(s,n)Li(n)dn, 


(3) 


where Li is the Legendre polynomial of order l. In practice 
we approximate eq. (3) by a discrete sum over the binned 
where we use As = 4 /i _1 Mpc and A/x=0.01 for every 
WiggleZ-BW and CMASS-BW region. We use the monopole 
(Z = 0) and quadrupole (Z = 2) of the 2-point functions, to 
analyse the redshift-space distortions, for separations s < 80 
h _1 Mpc. Our results are unchanged if large separations are 
used, whilst the increase in variance due to the finite number 
of mocks becomes significant. 

The covariance of each region is estimated from the mock 
WiZcola and BOSScola catalogs (see section 3.3). After cal¬ 
culating the covariances of the measurements in each overlap 
region from the COLA mock catalogs, we use inverse-variance 
weighting to obtain the ‘optimally combined’ measurements. 
For the statistic omb(s), the optimally combined function is 
calculated as 


N reg 

6,comb(^) — Ccomb ^ ^ Cj £z,z(s) (4) 

i=1 

where C CO mb is the overall covariance matrix, calculated from 
the estimations of the covariance matrices of individual regions 
Ci (see section 3.3). Results for the auto-2PCFs are shown in 
Figure 4, for individual regions (as lines) and for the combined 
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Figure 4. Combined monopole (top row) and quadrupole (bottom row) of £(s,/x) for the WiggleZ and CMASS auto-correlation function 
(left and middle columns) and the cross-correlation function (right column). Lines are measurements of individual regions, symbols display 
the combined measurements. Results are plotted as s 2 £z(s) as a function of separation s. 


measurements (as symbols). The different amplitude of clus¬ 
tering of the WiggleZ-BW and CMASS-BW galaxies reflects 
the difference in the type of halos these galaxies inhabit. Due 
to the limited volume where the correlations are measured, we 
correct our correlation function values by the ‘integral con¬ 
straint’ (Peebles 1980; Beutler et al. 2012). The corrections to 
the WiggleZ and BOSS correlations differ in each region and 
have values of the order of 8 x 10 -4 and 1 x 10 -3 respectively for 
the smaller regions (where the integral constraint is higher), 
and do not significantly affect the RSD model constraints. 


3.2 Cross-correlations between WiggleZ-BW and 
CMASS-BW clustering 


In addition to the auto-correlations, we also measured the 
cross-correlation between the two sets for tracers using the 
estimator 


Across (s, (l) 


DwDc — DwRc — RwDc + RwRc 
RwRc 


(5) 


where the W and C subscripts represent the quantities in 
the WiggleZ and CMASS galaxies, respectively. The cross¬ 
correlation function measurement provides an independent 
validation of the assumption that both galaxy types trace the 
same large structures on a range of scales, and also serves to 
test our linear and local galaxy bias model. 

To test the strength of the correlation between the tracers 
is, we also constrain the cross-correlation coefficient, r x , which 
is produced from the relation 


^WiggleZ x CMASS C 1 ®) ^X ( s )^Wigglez(s)^CMASs( s ) 5 (6) 






Figure 5. The cross-correlation coefficient r x for the WiggleZ- 
CMASS-BW correlations for each overlap regions, and when com¬ 
bining all regions, as a function of the smallest scale (in 

h -1 Mpc) of the fit. 


with |r x | ^ 1. On large scales in redshift-space, and assum¬ 
ing linear, deterministic bias, this quantity should tend to the 
value (Mountrichas et al. 2009): 

1 + | (Pw + Pc) + \PwPc ( V 

^ X, Kaiser — , — {•) 

y (l + §Pw + \Pw) (! + I Pc + \dc) 
Assuming bw = 1 and be = 2 (Reid et al. 2012; Contreras 
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Figure 6. Correlation matrix (normalised covariance matrix) of the 
WiggleZ-CMASS multipoles based on COLA mocks. The matrix is 
divided in the contributions from the monopole and quadrupole of 
WiggleZ-BW (W), CMASS-BW (C) and cross 2PCFS (X), and each 
pixel represents a separation s (in h~ 1 Mpc). 


et al. 2013), and a growth rate f(z — 0.54) = 0.75, when 
estimating f3 = f/b, it is expected that r x ,Kaiser = 0.997 ~ 1. 

We measure the value of r x from our data, assuming it is a 
constant on all scales (an assumption we do not expect to hold 
on scales smaller than 15 — 20 h~ x Mpc). Using the redshift 
space distortion model described in section 4.1, and the COLA 
mocks to build our covariance matrix, we use the correlation 
monopoles to fit for the bias parameters of the WiggleZ-BW 
and CMASS-BW galaxies and r x for each overlap region and 
the joint likelihood (see section 4.3 for details of the fitting 
procedure). 

Figure 5 presents the posterior probability distribution of 
r x , as a function of the minimum scale of fit Smin • Focusing 
on the fits to the combined regions, we can see that they are 
not consistent with r x = 1 at the 2cr level on scales Smin ~ 20 
hX 1 Mpc. This behaviour may be explained by a number of 
factors such as non-linear pairwise velocities, non-linear bias 
and stochasticity. CMASS galaxies tend to be hosted in the 
centres of large halos and in high density regions, precisely the 
regions that are avoided by WiggleZ galaxies. We expect that 
on large scales both galaxies trace similar structures, and this 
is confirmed in the measurements of r x being consistent with 
1 when fitting on large scales. 

Examining individual regions it can be noticed that it 
is region S22 which reduces the overall fit to r x . Its lower 
value of r x is driven by a high auto-correlation function in the 
WiggleZ-BW S22 region, although the scatter is compatible 
with the variance against mock catalogs. The best fits to the 
growth rate do not significantly change when the S22 region 
is excluded, and in the final fits we include all regions. 


3.3 Covariance estimation 

We estimate the correlations between the multipoles of the 
auto- and cross-2PCF by calculating the covariance matrix in 
each region n from COLA mocks. A deviation from the mean 
of a quantity X , in separation bin i, for the mock k can be 


written as 

At = X? - (X t ) (8) 

where, in our case, X corresponds to the monopole or 
quadrupole of the auto- or cross-2PCF in each bin. The co- 
variance matrix of each region n is determined as 

^mocks 

Cn,ij ~ T 7 (9) 

Smocks k~l 

After calculating C n for all regions, we can determine the com¬ 
bined covariance matrix (Blake et al. 2011a) 


= X 


( 10 ) 


Figure 6 shows the correlation matrix (normalised covariance 
matrix) for all our measurements, showing the strong correla¬ 
tion between the measurements of the two tracers. 

Since we used a large, but finite number of mock catalogs 
for the covariance estimation, there is an underestimation of 
the uncertainties. Following the work of Hartlap et al. (2007); 
Percival et al. (2014), we correct for the finite number of mocks 
by multiplying the variance estimated from the likelihood dis¬ 
tribution by 


1 + B(N hins - N p ) 
m<T 1 + 2A + N(N P + 1) 


( 11 ) 


where iVbins is the number of bins entering the fits, N p is the 
number of free parameters, and 

A — (-^mocks ^bins l)(^mocks ^bins 4), (f2) 

B — A(N moc ks — iVbms — 2), (13) 


where iV moc k s is the number of mock realisations. Also, the 
sample variance should be multiplied by 


X mocks 


mocks X\) ins 2 


(14) 


We use iVmocks = 480 and perform measurements in separation 
bins up to s = 80 h~ 1 Mpc. From constraining models using 
one-tracer auto-correlation function multipoles to simultane¬ 
ous fits using both auto and cross correlations, the m v factor 
lies in the range m v = 1.1 — 1.45. 


4 CONSTRAINTS ON COSMIC GROWTH 


4.1 Modelling the RSD 


Redshift-space distortions modify the 2-point clustering of 
galaxies on both large and small scales, which we will sum¬ 
marise here. Due to its peculiar velocity v, a galaxy at a po¬ 
sition in real space r gets mapped to s in redshift space: 


s = r + 


(1 + z)v ■ f „ 

—--- r 

H(z) 


(15) 


where r is the galaxy unit vector along the line of sight (LOS) 
direction, v r = v • r is the line-of-sight component of its veloc¬ 
ity, and H(z) is the Hubble parameter at a redshift z. 

On large scales, as described by Kaiser (1987), hereafter 
K87 (also see Hamilton 1998 for derivations in configuration 
space), matter overdensities 5m grow coherently as V • v oc 
— f5m where / = d\nG(a)/d\na is the linear growth rate 
of fluctuations. The evolution of the growth rate in certain 
models can be approximated by the evolution of the matter 
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density in the universe f(z) = ^m(z) 7 , where 7 = 0.55 in the 
case that the large-scale gravity obeys General Relativity; for 
alternative theories of gravity, 7 can take on different values 
(Linder & Cahn 2007). If we assume that the difference in 
clustering between dark matter and galaxies can be described 
by a linear bias model where S g — bSm then in Fourier space 
the redshift space galaxy overdensity takes the form 

<5 s (k) = (&+/7)Mk), (16) 

creating in configuration space the so-called ‘squashing’ effect 
on large scales. 

On small scales, in the non-linear regime for overdensi¬ 
ties and velocities, large structures appear elongated along 
the line of sight, creating the observed ‘Fingers of God’. In 
Fourier space this effect can be modelled by multiplying a 
Gaussian or a Lorentzian pairwise velocity distribution (i.e. a 
convolution of a Gaussian or exponential profile in configura¬ 
tion space) into the large-scale redshift-space distortion of the 
power spectrum. The simplest model, using Gaussian damping 
for the galaxy power spectrum in redshift space, is 

P s (k, p) = (b + f f i 2 ) 2 P m (k)e - (k ^ /H)2 (17) 

Pm(k) represents the non-linear real-space power spectrum 
and a v the pairwise velocity dispersion, which we approximate 
to be the same for all tracers is predicted to be 

a ' (z) = / Pee{k)dk (is) 

where, in the K87 formalism, Pee = Pm- However, this sim¬ 
ple model has been shown in simulations to be insufficiently 
accurate even on large scales, because there is not a perfect 
correlation between density and the velocity (divergence) field 
(e.g Okumura & Jing 2011; Kwan et al. 2012; de la Torre & 
Guzzo 2012; White et al. 2014). Scoccimarro (2004), hereafter 
S04, suggested a modification of the simple Kaiser formalism 
by including the velocity field terms. In the case of one tracer, 
the RSD in the galaxy auto- power spectrum reads: 

Pl(k,p) = [b 2 P S s(k) + 2 f i 2 fbPse(k) + f i 4 f 2 Pee(k)]e - (k ^ /H)2 

(19) 

where P 55 , Pse and Pee are the non-linear matter density- 
density, density-velocity and velocity-velocity power spectra, 
respectively. In our analysis, these terms are obtained from 
fitting formulae derived by Jennings (2012), from a suite of 
N-body simulations. In this case our fiducial model (based on 
WMAP5 results, see §1) predicts (via Eq.18) the large-scale 
velocity dispersion a v to be a v (z = 0.6) ~ 220 kms -1 . How¬ 
ever, we choose to leave a v as a free parameter to account 
for any additional non-linearities on smaller scales. Whilst 
there are additional improvements and implementations of 
RSD models (e.g. Taruya et al. 2010; Seljak & McDonald 2011; 
Reid & White 2011; Wang et al. 2013), this particular formal¬ 
ism has been successfully used in a number of studies (e.g. 
Blake et al. 2011b; de la Torre et al. 2013; Blake et al. 2013), 
and, as we will see below, reproduces the expected constraints 
on the growth rate from the COLA mock catalogs and pro¬ 
vides a good description of the galaxy anisotropic clustering 
at the current statistical level. 

In the case of the redshift-space cross-power spectrum, as¬ 
suming that both tracers are described by the same dispersion 
parameter a v , we can write the large-scale terms as 

Pl(k, At) = [bib 2 Pss(k) + p 2 f(bi + b 2 )Pse{k) + 

p A f 2 Pee{k)\ x e ~( k ^v/H(z)) 2 (20) 


where b\ and b 2 are the biases of the different tracers. Since we 
are measuring the multipoles of the 2PCF, we calculate first 
the power spectrum multipoles as 

PHk) = ^+1 f 1 P’(k,riL l (jI)dp, ( 21 ) 

where l is the multipole order and Li is the Legendre polyno¬ 
mial of order l. Then for the 2-point correlation function in 
configuration space we have 

tl (s) = ^ J Pi (k)j i (ks)k 2 dk ( 22 ) 

where ji is the spherical Bessel function of order l. 


4.2 Tests using COLA mocks 

We tested the validity of these models using our COLA mock 
catalogs. In summary, we compared the K87 and S04 mod¬ 
els for the large-scale distortions to P(k) (calculated using 
our fiducial cosmological parameters), using a Gaussian func¬ 
tion for the small-scale damping (we also tried fits using the 
Lorentzian profile without significant differences), constrain¬ 
ing the growth rate / at the simulation output redshift z = 0.6, 
marginalising over the bias of each tracer and the common ve¬ 
locity dispersion a v . We performed these fits for every COLA 
mock on scales s < 80 /i -1 Mpc, although changes when using 
larger scales were not significant. 

Figure 7 shows the mean and standard deviation of the 
best-fitting values of f(z — 0.6) across the mocks. For the 
WiZcola mocks the K87+Gaussian model tends to underpre¬ 
dict the value of the growth rate whereas the S04+Gaussian 
model agrees well for scales s m i n > 20 /i -1 Mpc. For the BOSS- 
cola mocks the differences are less pronounced, but the input 
growth rate is recovered with a systematic error less than the 
statistical error. In both cases the goodness of fit is similar 
with x 2 /d.o.f.~ 1 for both WiggleZ and BOSS COLA mocks 
on larger scales s m in > 20 /i -1 Mpc, worsening considerably 
on scales s m in < 10 /i _1 Mpc. In what follows we will use the 
S04+Gaussian model for our parameter fits. 

There are, however, specific differences in the scale of va¬ 
lidity of the models depending on which tracer is used. It can 
be seen that for low-bias galaxies represented by the WiZ¬ 
cola mocks, the agreement between the model fits and the 
input value of f(z = 0.6) extends to lower scales than in the 
case of galaxies residing in more massive halos. Although the 
Kaiser effect is stronger for lower bias galaxies, the higher 
non-linearities arising from the formation and high-clustering 
of high-mass halos lead to a model break-down on larger 
scales. For the particular case of the WiggleZ-BOSS overlap, 
Smin — 24 h ~ x Mpc is the minimum scale where both mod¬ 
els recover adequately the fiducial growth rate with negligible 
systematic error. 

Multitracer approach 

Having chosen the model to analyse redshift-space cluster¬ 
ing, we examined the consequences of using multiple tracers 
when recovering model parameters. For each realisation of the 
COLA mocks we fit the S04+Gaussian RSD model first using 
the autocorrelations independently, then analysing both auto¬ 
correlations but considering the common covariance matrix, 
and lastly adding the cross-correlations using the monopole 
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Figure 7. Constraints on the linear growth rate f(z = 0.6) us¬ 
ing COLA mock galaxies with the same angular and redshift selec¬ 
tion functions as the galaxy survey data. Fits to the growth rate 
f(z = 0.6) (= 0.76, for the WMAP5 cosmology, green solid line) are 
performed for the models as a function of Smim with Smax = 80 
h -1 Mpc. Colored shades cover the lcr (68%) confidence interval for 
/, defined using the dispersion in the best-fitting values across the 
480 mocks. 


and quadrupole of £(s,/i) in the range 24 < s < 80 /i _1 Mpc. 
Results are shown in Figure 8, which displays the lcr contours 
enclosing the best-fitting values of f(z = 0.6) and a v . The 
expected values for these parameters are consistent with our 
COLA constraints, and the approximation that a v is the same 
for both tracers is valid for this range of scales. The constraints 
for the parameters are correlated between the two surveys, 
with a cross-correlation coefficient p^ wc — cr wc/ <7 w cr c ~ 0-4 
for the growth rate between both surveys. 

When analysing the 2PCFs of the two tracers simultane¬ 
ously, taking into account the common covariance, an improve¬ 
ment in the measurement of / is obtained, of the order 30% 
compared to using the BOSScola mocks alone (which because 
of a higher bias, have a lower value of [3 and hence a lower 
signal) but only 5% compared to using WiZcola mocks alone. 
Adding the cross-2PCF produces an improvement of 20% com¬ 
pared to the WiZcola-only constraints, mostly due to an in¬ 
creased signal in the shot-noise dominated regime. Analysing 
individual mocks shows that the improvement also varies in 
each realisation. As predicted by McDonald X Seljak (2009), 
Gil-Marm et al. (2010) and Blake et al. (2013), although our 
tracers have big differences in their biasing, due to the sparsity 
of our sampling we are in the regime where shot noise domi¬ 
nates and improvement via the cancellation of cosmic variance 
is small. 


4.3 Data Fitting procedure 

In our analysis we fixed the cosmological parameters of the 
matter power spectra to the best-fit WMAP5 model (Komatsu 
et al. 2009), the fiducial cosmology of our COLA mocks, and 
constrain the parameters (bw,bc,f(z = 0.54), a v ). Due to 



f(z=0.6) 


Figure 8. Constraints for RSD parameters in the S04 model for 
fits to COLA mocks in the range 24 < s < 80 /i -1 Mpc using 2PCF 
multipoles. Results for the full fit and for individual subsamples, 
with the contours enclosing 68% of the best-fitting parameters are 
shown. The four cases considered are WiZcola mocks alone (blue 
solid line), BOSScola mocks alone (red dashed line), joint analysis 
of the auto-correlations (black solid line), and an analysis further 
adding the cross-correlation function (green dashed line). 


the degeneracy of the first three parameters with cr 8 (z), the 
r.m.s. of the matter density field in 8 /i _1 Mpc spheres, we are 
effectively constraining (bwos^bcos, fcrg, <J V ). When we also 
include the WiggleZ-CMASS cross-correlation in the analysis, 
we additionally fit for the parameter r x • We compare the con¬ 
straints from the single-tracer model for each galaxy type to 
each other, and then include the common covariance and the 
cross-correlations in the cosmological fits. 

We use the monopole and quadrupole of the tracers, and 
present results as a function of the minimum-scale fitted s m i n , 
with Smax = 80 /i _1 Mpc. We execute a Maximum Likelihood 
parameter estimation test, where we minimise the quantity 


i 0 

(23) 

where X is one of the elements of the vector formed by the 
multipoles of £(s,/x) of WiggleZ-BW, CMASS-BW and/or 
WiggleZ-CMASS-BW correlations. We explore the parameter 
space using a Monte Carlo Markov Chain method (MCMC) 
imposing the prior that all parameter values must be bigger 
than zero. 


4.4 Fits for the growth rate 

Figure 9 presents the parameter fits of /cr 8 (z = 0.54) fitting 
the monopole and quadrupole of the WiggleZ and CMASS 
auto- and cross-correlation on scales between Smin — 24 
h -1 Mpc, and s macc = 80 h~ x Mpc. As shown in the previ¬ 
ous section, s m in = 24 /i _1 Mpc is the minimum scale where 
there are not important systematic deviations in the param¬ 
eters from the study with the COLA mocks, and our fits to 
the observed data follow this trend. Table 2 lists the results 
for the parameter fits. 

Comparing the single-tracer fits for WiggleZ-BW and 
CMASS-BW galaxies, there is agreement at the lcr level for 
the values of fcrs(z — 0.54), meaning that when fitting to 
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Figure 9. Fits to the RSD model parameters using correlation mul¬ 
tipoles with 24 < s < 80 h -1 Mpc. We show results when analysing 
individual surveys and joint constraints. The purple line shows the 
prediction from WMAP5 cosmology. 


Figure 10. Fits to the growth rate fas(z) from different galaxy 
surveys: 6dF (Beutler et al. 2012), 2dFGRS (Hawkins et al. 2003), 
GAMA (Blake et al. 2013), WiggleZ (Blake et al. 2011b), SDSS 
LRGs (Samushia et al. 2012), CMASS-DR9 (Reid et al. 2012), 
VVDS(Guzzo et al. 2008), and VIPERS (de la Torre et al. 2013). 


these scales there is evidence of no systematics depending on 
the type of galaxy used. Our constraints on the growth rate are 
consistent with our fiducial cosmology fas(z m 0.54) = 0.46. 

Consistent with previous work, we recover that the bias 
of the WiggleZ-BW galaxies, bw ~ 1, is smaller than that 
of the CMASS-BW galaxies, be ~ 2. The value of the best¬ 
fitting chi-squared statistic indicates that the model provides 
a reasonable fit to the data in all cases. For the pairwise dis¬ 
persion, values for the different tracers are consistent with the 
predicted value from theory (section 4.1). 

Combining the two tracers including their cross¬ 
covariance yields slightly better constraints for fas(z = 0.54) 
at the level of 10% (compared to WiggleZ constraints alone). 
This result indicates that for these tracers, in a low den¬ 
sity regime (where the common cosmic variance cancellation 
does not improve the constraints, see Blake et. al 2013), even 
in the presence of a slightly larger Hartlap-Percival correc¬ 
tion, the improvement is due to reduced shot noise. When 
including the cross-correlations the improvement is of the or¬ 
der of 20% (again, compared with WiggleZ constraints alone). 
In the case when we include the cross-correlations, we ob¬ 
tain our poorest value for x 2 /d.o.f., implying that our sim¬ 
ple constant r x model may not describe all of the complex¬ 
ities of the cross-correlation. Given this result, we quote as 
result of our paper for the growth rate constraint the one ob¬ 
tained when we combine only the auto-correlations, yielding 
f(7s{z = 0.54) = 0.413 zb 0.054. 


5 SUMMARY & CONCLUSIONS 

In this work we have presented the first cosmological RSD 
analysis using data from two overlapping surveys, WiggleZ 
and CMASS. After defining the overlap volumes, we mea¬ 
sured 2-point auto- and cross-correlations functions of these 
tracers; after obtaining their multipoles and calculating their 
cross-covariance using N-body mock catalogs, we compared 


them with RSD models in order to measure the growth rate 
of structure fas at an effective redshift z — 0.54. Our main 
findings are: 

• The cross-correlation coefficient r x between the WiggleZ- 
BW and CMASS-BW galaxies agrees with the expectation 
that on large scales, the two classes trace similarly the large 
scale structure with r x ~ 1. On smaller scales s < 20 /i _1 Mpc, 
r x < 1, likely produced by a combination of a number of 
factors such as non-linear pairwise velocities, non-linear bias 
and stochasticity. 

• We tested redshift-space distortion models in mock cata¬ 
logues simulating WiggleZ and CMASS galaxies, including the 
selection functions of our overlapping volumes. When fitting 
scales s > 24 /i _1 Mpc we recover our fiducial cosmological 
parameters using different tracers, and that a single velocity 
dispersion provides an adequate description for the distortions 
in our range of scales. We confirmed a lack of a significant im¬ 
provement when using the multitracer technique, given the 
sparsity of the sampling for these tracers. 

• The fits to fag (z) from all tracers are consistent with each 
other and with the predictions of a ACDM universe, showing 
no evidence for strong modelling systematic errors as a func¬ 
tion of galaxy type. 

As shown in Figure 10, our combined fit for the growth 
rate fag(z = 0.54) = 0.413 ± 0.054 is in excellent agreement 
with estimates from different surveys. Although more sophis¬ 
ticated models for the RSD can be employed, the motivation 
for our work was to show consistency in the cosmological fits 
when using different tracers. This agreement provides further 
strong evidence for the robustness in the growth rate mea¬ 
surements which are important for answering the outstanding 
questions on the nature of dark energy and large-scale gravity. 
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Table 2. Fits to the RSD model parameters when using multipoles with 24 < s < 80 h 1 Mpc assuming WMAP5 <js(z = 0) = 0.812 


Tracers 

b\yri8(z = 0.54) 

b c o- 8 (z = 0.54) 

MsCz = 0.54) 

a v (kms 1 ) 

r X 

X 2 /d.o.f 

d.o.f. 

WiggleZ only 

0.651T0.046 

- 

0.409±0.059 

205±144 

- 

1.11 

28 - 4 

CMASS only 

- 

1.204T0.062 

0.466±0.074 

130±116 

- 

1.43 

28 - 4 

WiggleZ+CMASS 

0.646T0.043 

1.233 ±0.054 

0.413±0.054 

117±113 

- 

1.28 

56 - 5 

WiggleZ±CMASS±x2PCF 

0.648T0.038 

1.242±0.043 

0.403±0.048 

88±104 

0.93 ± 0.03 

1.57 

84 - 6 
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